Correlations between kHz QPO and Low Frequency Features 
Attributed to Radial Oscillations and Diffusive Propagation in the 
Viscous Boundary Layer Around a Neutron Star 

Lev Titarchuk 

NASA/ Goddard Space Flight Center, Greenbelt MD 20771, and George Mason 
University /CSI, USA; titarchuk@lheavx.gsfc.nasa.gov 

Vladimir Osherovich 
NASA/Goddard Space Flight Center/RITSS, Greenbelt MD 20771 USA; 

vladimir @urap . gsf c . nasa. gov 

Received ; accepted 



- 2 - 
ABSTRACT 

We present a dimensional analysis of two characteristic time scales in the 
boundary layer where the disk adjusts to the rotating neutron star (NS). The 
boundary layer is treated as a transition region between the NS surface and the 
first Keplerian orbit. The radial transport of the angular momentum in this 
layer is controlled by a viscous force defined by the Reynolds number, which 
in turn is related to the mass accretion rate. We show that the observed low- 
Lorentzian frequency is associated with radial oscillations in the boundary layer, 
where the observed break frequency is determined by the characteristic diffusion 
time of the inward motion of the matter in the accretion flow. Predictions of our 
model regarding relations between those two frequencies and frequencies of kHz 
QPO's compare favorably with recent observations for the source 4U 1728-34. 
This Letter contains a theoretical classification of kHz QPO's in NS binaries 
and the related low frequency features. Thus, results concerning the relationship 
of the low-Lorentzian frequency of viscous oscillations and the break frequency 
are presented in the framework of our model of kHz QPO's viewed as Keplerian 
oscillations in a rotating frame of reference. 

Subject headings: accretion, accretion disks — diffusion — stars:individual (4U 
1728-34, Sco X-1) — stars:neutron — X-ray:star — waves 
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1. Introduction 

The discovery of kilohertz quasiperiodic oscillations (QPO's) in the low mass X-ray 
neutron star (NS) binaries (Strohmayer et al. 1996; Van der Klis et al. 1996 and Zhang et 
al. 1996) has stimulated both theoretical and observational studies of these sources. In the 
upper part of the spectrum (400- 1200 Hz) for most of these sources, two frequencies and 
Vh have been seen. Initially, the fact that for some sources, the peak separation frequency 
Az/ = Vh — i^k docs not change much led to the beat frequency interpretation (Strohmayer et 
al. 1996; Van der Klis 1998) which was presented as a concept for the first time in the paper 
by Alpar & Shaham (1985). Beat-frequency models, where the peak separation is identified 
with the NS spin rate have been challenged by observations: for Sco X-1, Ai/ varies by 40% 
(van der Klis et al. 1997 hereafter VK97) and for source 4U 1608-52, Ai/ varies by 26% 
(Mendez et al. 1998). Mounting observational evidence that Ai/ is not constant demands a 
new theoretical approach. For Sco X-1, in the lower part of the spectrum, VK97 identified 
two branches (presumably the first and second harmonics) with frequencies 45 and 90 Hz 
which slowly increase in frequency when and Uh increase. Furthermore, in the spectra 
observed by Rossi X-ray Timing Explorer (RXTE) for 4U 1728-34, Ford and van der Khs 
(1998, herein FV98) found low frequency Lorentzian (LFL) oscillations with frequencies 
between 10 and 50 Hz. These frequencies as well as break frequency, v^yreak of the power 
spectrum density (PSD) for the same source were shown to be correlated with and Vh. It 
is clear that the low and high parts of the PSD of the kHz QPO sources should be related 
within the framework of the same theory. Difficulties which the beat frequency model faces 
are amplified by the requirement of relating the observed low frequency features, described 
above, with and Vh. 

Recently, a different approach to this problem has been suggested: kHz QPO's in the 
NS binaries have been modeled by Osherovich & Titarchuk (1999) as Keplerian oscillations 
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in a rotating frame of reference. In this new model the fundamental frequency is the 
Keplerian frequency (the lower frequency of two kHz QPO's) 



where G is the gravitational constant, M is the NS mass, and R is the radius of the 
corresponding Keplerian orbit. The high QPO frequency Uh is interpreted as the upper 
hybrid frequency of the Keplerian oscillator under the influence of the Coriolis force 



where Q is the angular rotational frequency of the NS magnetosphere. 

For three sources (Sco X-1, 4U 1608-52 and 4U 1702-429), we demonstrated that the 
solid body rotation (Q = Qq = const) is a good first order approximation. Slow variation 
of O as a function of i^k within the second order approximation is related to the differential 
rotation of the magnetosphere controlled by a frozen-in magnetic structure. This model 
allows us to address the relation between the high and low frequency features in the 
PSD of the neutron systems. We interpreted the ~ 45 and 90 Hz oscillations as 1st and 
2nd harmonics of the lower branch of the Keplerian oscillations in the rotating frame of 
reference: 



where 5 is the angle between Q and the vector normal to the plane of the Keplerian 
oscillations. For Sco X-1, we found that the angle S — 5.5° fits the observations. 

In this Letter we include the LFL oscillations and related break frequency phenomenon 
in our classification. We attribute LFL oscillations to radial oscillations in the viscous 
boundary layer surrounding a neutron star. According to the model of Shakura & Sunyaev 
(1973, hereafter SS73), the innermost part of the Keplerian disk adjusts itself to the 
rotating central object (i.e. neutron star). The recent modelling by Titarchuk, Lapidus & 
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i/L = {fl/n){uk/iyh) sin (5, 
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Muslimov (1998, hereafter TLM) led to the determination of the characteristic thickness 
of the viscous boundary layer L. In the following section, we present the extension of this 
work to relate the frequency of the viscous oscillations and inbreak with i/^. Comparison 
with the observations is carried out for 4U 1728-34. The last section of this Letter contains 
our theoretical classification of kHz QPO's and related low frequency phenomena. 



2. Radial Oscillations and Diffusion in the Viscous Boundetry Layer 

We define the boundary layer as a transition region confined between the NS surface 
and the first Keplerian orbit. The radial motion in the disk is controlled by the friction 
and the angular momentum exchange between adjacent layers resulting in the loss of the 
initial angular momentum by an accreting matter. The corresponding radial transport of 
the angular momentum in a disk is described by the equation (e.g. SS73): 

^Jr^""^'^ = 

where M is the accretion rate, and Wrip is the component of a viscous stress tensor which is 
related to the gradient of the rotational frequency uj, namely 

Wr^ = -'^VHR^, (5) 

where H is a half-thickness of a disk, and r) is the turbulent viscosity. The nondimensional 
parameter which is essential for equation (4) is the Reynolds number for the accretion fiow 

_ M _ 3RVr 

which is the inverse a— parameter in the SS73-model; Vr is a characteristic velocity, Vt and 
If are a turbulent velocity and related turbulent scale respectively. 

Equations cu — cuq at R = Rq (NS radius) and a; = a;K at R — Rout (radius where the 
boundary layer adjusts to the Keplerian motion), and ^ = ^ at R = Rout were assumed 
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by TLM as boundary conditions. Thus the profile u{R) and the outer radius of the viscous 
boundary layer R^ut are uniquely determined by these boundary conditions. 

Presenting u!{R) in terms of dimensionless variables: namely angular velocity 9 — cu/cuq, 
radius r = R/Rq {Rq = xqRs, Rg = 2GM/c^ is the Schwarzschild radius), and mass 
m = M/Mq, we express Keplerian angular velocity as 

9k = G/iaKv'/'), (7) 

where Uk = m(a;o/3)'^''^(z/o/363 Hz) and the NS rotational frequency Uq has a particular 
value for each star. The particular coefficient, 6, presented in formula (7) is obtained for 
the frequency of nearly coherent (burst) oscillations for 4U 1728-34, i.e. for uq — 363 Hz. 
The solution of equations (4-5 ) satisfying the above boundary conditions is 

9{r) = L>ir-^ + (1 - L>i)r-', (8) 

where Di = {9out - r^t)/{r'^t - r^t) and 9out = 9K{rout)- Equation 9' (rout) = 9'K{rout) 
determines Tout- 

\eout - ^i7r J + 2(1 - D,)r^,. (9) 

The solution of equations (4-5) subject to the inner sub-Keplerian boundary condition 
has a regime corresponding to the super-Keplerian rotation (TLM). For such a regime 
matter piles up in the vertical direction thus disturbing the hydrostatic equilibrium. The 
vertical component of the gravitational force prevents this matter from further accumulation 
in a vertical direction and drives relaxation oscillations. The radiation drag force, which is 
proportional to the vertical velocity, determines the characteristic decay time of the vertical 
oscillations (TLM). The characteristic time t^i over which the matter moves inward through 
this region, bounded between the innermost disk and relaxation oscillations zone is 

tr ~ -, (10) 
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where L — — Rq is the characteristic thickness of this region. Even though the specific 
mechanism providing the modulation of the observed X-ray fiux over this timescale needs 
to be understood, this timescale apparently "controls" the supply of accreting matter into 
the innermost region of the accretion disk. Any local perturbation in the transition region 
would propagate diffusively outward over a timescale 



where Ifp is the mean free path of a particle. 

Note, that the 7— parameter is proportional to the accretion rate (see Eq. 6), and 
therefore Vr oc 7. Using this relationship, we can exclude Vr from the above equations and 
get the relations for the corresponding inverse timescales (frequencies). For the frequency 
of viscous oscillations 

7 



font ~ fO 

and for the break frequency, related to the diffusion 



(12) 



In the following section, we compare the predictions of this model with the observations 
and also establish the theoretical relation between i/y and lybreak- 



3. Compcirisons with Observations 

The results of FV98 for the low frequency Lorentzian in the X-ray binary 4U 1728-34 
are presented in Figure 1 and for the break frequency lybreak in Figure 2. In Figure 1, crosses 
represent the frequencies (with the appropriate error bars) observed during four days. Data 
collected on February 16 (open circles) are situated apart from the rest of observations and 
they are not included in the empirical power law fit which is suggested by FV98. In the work 
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discussed above, the authors plotted the observed low frequencies versus high-frequency 
QPO which for all days, except February 16 was Vk and apparently for February 16 it was 
Uh. Our theoretical curve for versus Uk is based on equation (12). The dependence on 
this parameter is rather strong: the parabola — 38024 — 73076 ■ ak + 35732 • a| has a 
minimum at = 1.03, which determines the best fit. Using Q/27r = 340 Hz in the upper 
hybrid relation (2), we calculate i^k for the points observed on February 16 and show that 
they belong to the set of frequencies modeled by our theoretical curve for the viscous radial 
oscillations (closed circles). 

Identification of the observed Vbreak with the inverse diffusion time (formulas 11 and 13) 
is illustrated by theoretical curves in Figure 2. It is worth noting that these two correlations 
with kHz frequencies are fit by two theoretical curves using only one parameter a^. The 
X^— dependence on is obtained with inclusion of all data points for the break and low 
frequency correlations (75 data points). The theoretical dependences of Uk and rout on 
7— parameter are calculated numerically using equations (7) and (9) and employed here for 
calculations of the theoretical curves in Figures 2 and 3 using equations (12) and (13). 

We were unable to interpret data for Vbreak collected in February 16 (open circles). 
Neither z/„ nor Uhreak in our theory have a power law relation with z/^. However, the 
theoretical relation between lybreak and i/y, shown in Figure 3 by a solid curve, is close to the 
straight line (in log-log diagram), suggesting an approximate power law 

Ubreak^OMlu'/K (14) 

This relation is derived from the theoretical dependence for the best fit parameter = 1.03. 
Observations of FV98 (except February 16) are also presented in the Figure 3. 

4. Discussion and Conclusions 

We present a model for the radial oscillations and diffusion in the viscous boundary layer 
surrounding the neutron star. Our dimensional analysis has identified the corresponding 
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frequencies and inbreak which are consistent with the low Lorentzian and break frequencies 
for 4U 1728-34. and predicted values for Uhreak related to the diffusion in the boundary layer 
are consistent with the break frequency observed for the same source. Both oscillations 
(Keplerian and radial) and diffusion in the viscous boundary layer are controlled by the 
same parameter - Reynolds number 7 which in turn is related to the accretion rate. It is 
shown in TLM that i^k is a monotonic function of 7. Therefore, the observed range of v^, 
(350-900 Hz) corresponds to the range 1 < 7 < 5 (or 0.2 < o; < 1). 

The results in this Letter extend the classification of kHz QPO's and the related low 
frequency phenomena suggested by Osherovich & Titarchuk 1999. Figure 4 summarizes the 
new classification. Solid lines represent our theoretical curves and open circles observations 
for Sco X-1 (from VK97). As one can see, formulas (2) and (3), for the Keplerian 
oscillator under the infiuence of the Coriolis force, reproduce the observations well. Indeed, 
Au — i/fi — Uk is not constant, as observed (see OT99 for details of comparisons of the 
data with the theory). Effectively, the main viscous frequency and the diffusive inbreak 
introduce the second oscillator with two new branches in the lower part of the spectra. The 
unifying characteristic of spectra for both oscillators is the strong dependence on 1/^. This 
common dependence on Uk can be viewed as a result of the interaction between Keplerian 
oscillator and the viscous oscillator which share the common boundary at the outer edge of 
the viscous transition layer. 

Our parametric study indicates that the power law index 1.6 in Eq. (14) should be the 
same for different neutron stars. We expect a similar relation for black holes but with a 
distinctly different index. 

The found value of leads ultimately to independent constraints in the determination 
of mass and radius for the neutron star (Haberl & Titarchuk 1995). 
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Fig. 1. — Main Viscous Frequency (crosses) versus the frequency of the kHz QPOs for 4U 
1728-34 (FV98). Open circles are for the data of February 16. Sohd circles are for the data of 
16th of February recalculated using the Keplerian frequencies with the rotational frequency 
340 Hz (see text). The sohd line is the best theoretical fit for the main viscous frequency 
calculated using Eqs. (4-5, 9, 12) corresponding to the parameter = 1.03. 



Fig. 2. — Break Frequency of a broken power law function of the power spectrum density 
(crosses) versus the frequency of the kHz QPOs for 4U 1728-34 (FV98). Open circles are 
for the data of February 16. The solid line is the theoretical curve for break frequencies 
calculated using Eqs. (4-5, 9, 13) corresponding to the parameter = 1.03. 



Fig. 3. — Break frequency versus main viscous frequency (stars) for all data points except 
those of February 16. Solid line is the theoretical curve which is almost a power law with 
index 1.61. 



Fig. 4. — Classification of the QPO frequencies in kHz QPO sources. Open circles of QPO 
data by VK97. 
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